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Abstract 

The mechanism regulating expression of late genes, encoding viral structural components, is an unresolved problem in the 
biology of DNA tumor viruses. Here we show that BGLF4, the only protein kinase encoded by Epstein-Barr virus (EBV), 
controls expression of late genes independent of its effect on viral DNA replication. Ectopic expression of BGLF4 in cells 
lacking the kinase gene stimulated the transcript levels of six late genes by 8- to 10-fold. Introduction of a BGLF4 mutant 
that eliminated its kinase activity did not stimulate late gene expression. In cells infected with wild-type EBV, siRNA to BGLF4 
(siG4) markedly reduced late gene expression without compromising viral DNA replication. Synthesis of late products was 
restored upon expression of a form of BGLF4 resistant to the siRNA. Studying the EBV transcriptome using mRNA-seq during 
the late phase of the lytic cycle in the absence and presence of siG4 showed that BGLF4 controlled expression of 31 late 
genes. Analysis of the EBV transcriptome identified BGLF3 as a gene whose expression was reduced as a result of silencing 
BGLF4. Knockdown of BGLF3 markedly reduced late gene expression but had no effect on viral DNA replication or 
expression of BGLF4. Our findings reveal the presence of a late control locus encompassing BGLF3 and BGLF4 in the EBV 
genome, and provide evidence for the importance of both proteins in post-replication events that are necessary for 
expression of late genes. 
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introduction 

Late genes encode structural proteins necessary for virion 
assembly. A common theme among DNA viruses is the strict 
dependence of late gene expression on the onset of viral DNA 
replication. Disruption of replication, using inhibitors or mutating 
a replication-essential gene, blocks synthesis of late products. Tlie 
link between these two processes led to models that focus on 
genome amplification as the principal regulator of late gene 
expression. These models propose that changes in DNA modifi- 
cations such as a decrease in methylation as a result of de novo 
DNA synthesis, or displacement of viral or cellular repressors 
bound to elements in late promoters by the replication machinery 
trigger late gene expression. However, mechanisms that link late 
gene expression to replication have not been elucidated. Reports 
in herpesviruses suggest that replication per se is not sufficient to 
activate late gene expression. Pioneering work by Ren Sun and 
colleagues in murine herpesvrrus-68 (MHV-68) identified four 
early viral proteins, ORF18, 24, 30 and 34, to be required for 
expression of late genes but dispensable for viral DNA replication 
[1,2,3,4]. Homologs of ORFs 18, 24 and 34 in human 



cytomegalovirus (hCMV) map to three unique long (UL) 
sequences 79, 87 and 95, respectively [5,6]. The function of the 
MHV68 and hCMV proteins in activating late gene expression 
has not been elucidated. In Epstein-Barr virus (EBV), the only 
protein so far characterized as essential for activation of late genes 
and not DNA replication is BcRFl, a homolog of ORF24 in 
MHV-68 and UL87 in CMV [7]. BcRFl is a TATA box binding- 
like protein that specifically binds to a non-canonical TATA 
element (TATT) present in most late promoters [8,9] . 

Viral late promoters differ from viral promoters of other kinetic 
classes and cellular promoters that rely on transcription factor- 
binding sites located upstream of the TATA box. Activation of late 
promoters is primarily dependent on a distinct TATA box and a 
downstream initiator sequence that spans the transcription start 
site (TSS) [9,10,11]. Involvement of upstream elements in 
regulation of late promoters is postulated to modulate transcription 
efficiency [12,13,14,15]. 

BGLF4 is the only Ser/Thr protein kinase encoded by EBV 
[16]. Homologs of BGLF4 are conserved among the herpesvirus 
family [17]. These conserved human herpesvirus protein kinases 
(CHPKs) share sequence and positional similarities, but exhibit 
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Author Summary 

Epstein-Barr virus (EBV) is linked to the development of 
several types of cancer. Synthesis of structural proteins, a 
group of proteins that forms the protein shell around the 
viral genome, is essential for EBV infection and pathogen- 
esis. Genes encoding structural proteins are expressed late 
in the viral life cycle after amplification of the viral genome. 
The mechanism controlling expression of this group of 
proteins represents a longstanding conundrum in EBV and 
other DNA viruses. In this report, we demonstrate that two 
EBV regulatory proteins control synthesis of mRNAs 
encoding viral structural proteins. These two proteins 
are: BGLF4, a protein kinase conserved in all herpesviruses, 
and BGLF3, a protein of unknown function with no cellular 
counterparts. We present evidence that the enzymatic 
activity of BGLF4 is required after replication of viral DNA 
to stimulate expression of structural proteins. BGLF3 and 
BGLF4 are expressed from the same locus in the genome; 
the two proteins work in concert and independently to 
promote expression of viral genes encoding structural 
proteins. Our findings provide novel insights into control 
of expression of genes encoding viral structural proteins. 
The enzymatic activity of BGLF4 is a potential target for 
development of new antiviral drugs. 



both unique and overlapping substrate specificities [18,19]. 
BGLF4, like other CHPKs, functionally mimics cellular cyclin- 
dependent kinases but displays broader substrate specificity 
[17,19,20,21,22]. For instance, protein array phosphorylation 
experiments identified 21 EBV proteins as putative substrates of 
BGLF4; half of these proteins are shared targets with CDKl/ 
cyclinB [23]. A number of cellular cyclin-dependent kinase 
substrates are also modified by BGLF4 such as pRB, p27, 
condensin, MCM4, stathmin, elongation factor 1 delta, and 
nuclear lamin A/C [19,20,21,22,24,25,26,27,28,29]. Phosphory- 
lation of lamin A/C, another CDK substrate [30], by BGLF4 
causes dissolution of the nuclear lamina, a step considered essential 
for egress of viral capsids from the nucleus [27,31,32]. 

BGLF4 has been extensively investigated to determine its 
potential role in viral DNA replication. Expression of BGLF4 
occurs during the early phase of the lytic cycle and reaches maximal 
levels following viral DNA replication [33]. BGLF4 localizes to 
replication compartments, a site of viral DNA replication and late 
gene expression [34,35]. BGLF4 modulates gene expression by its 
capacity to phosphorylate transcription factors, histones and 
chromatin modifying enzymes, such as BMRFl, ZEBRA, EBNA2, 
EBNA-LP, IRF3, UXT, HDACl, HI, and TIP-60 
[16,36,37,38,39,40,41,42,43,44,45]. However, the evidence in favor 
of a role of BGLF4 in regulating early gene expression and viral 
DNA replication during the EBV lytic cycle is conflicting. 
Abolishing expression of BGLF4 had no significant effect on early 
gene expression [46] and had littie (<2-fold) to no impact on viral 
DNA replication [28,31,32,44,45,47]. Perhaps the most established 
phenotype for disrupting the expression or the activity of BGLF4 is 
the significant reduction in the amount of virus released by lytic 
infected cells [27,31,48]. This phenotype is consistent with a crucial 
role for BGLF4 in events occurring after viral DNA replication. 

In this study, we employed knockout and knockdown 
approaches to thoroughly investigate the role of BGLF4 in 
regulating lytic viral gene expression and its correlation with DNA 
replication. We found that the kinase activity of BGLF4 is 
necessary for expression of late genes but not viral DNA 
replication. Using RNA-seq analysis we demonstrated that 



expression of most late genes is reduced in the absence of BGLF4. 
We identified a late control locus encompassing BGLF4 and 
BGLF3. Knockdown of BGLF3 abolished late gene expression 
independent of any effect on the level of BGLF4 or viral DNA 
replication. Our findings identify two EBV genes as regulators of 
late gene expression and establish the presence of additional 
checkpoints beyond viral DNA replication that are necessary for 
expression of EBV structural proteins. 

Results 

BGLF4 significantly up-regulates expression of the late 
BFRF3 protein 

BGLF4 is transcribed as a multicistronic message that also 
encodes BGLF5, a DNA alkaline exonuclease involved in host cell 
shutoff [47,49]. To study the role of BGLF4 in regialating the lytic 
cascade we used delta G4/ G5 cells, which are 293 cells harboring 
BGLF4/BGLF5-null EBV bacmid virus [47]. Cells were harvested 
48 h after transfection to allow for analysis of viral DNA 
replication and late gene expression. Typically, transient expres- 
sion of ZEBRA in EBV positive cells induces expression of three 
major kinetic classes of lytic genes, these are: very early, early and 
late genes [50,5 1] . In delta G4/G5 cells, we found that ZEBRA 
activated the EBV lytic cycle as demonstrated by expression of the 
early bmrfl gene encoding the DNA polymerase processivity 
factor (Fig, lA lane 2). However, expression of the late BFRF3 
protein (FR3), the minor viral capsid protein, was noticeably low, 
only 2.5-fold higher in cells transfected with ZEBRA relative to 
empty vector (CMV) (Fig. lA, lanes 1 and 2). Complementation 
with BGLF4 significantly enhanced the level of FR3 to 20.5-fold 
relative to empty vector and 8.2-fold relative to ZEBRA- 
transfected cells (Fig. lA, lane 3). Complementation with BGLF5 
did not enhance expression of FR3, a result suggesting that the 
effect on late gene expression is specific to BGLF4. BMRFl is a 
bona fide substrate of the BGLF4 kinase [16,36,52]. Therefore, 
hyper-phosphorylation of BMRFl is a marker for the intracellular 
kinase activity of BGLF4. Hyper-phosphorylated BMRFl was 
only detected in lytic cells transfected with the BGLF4 expression 
vector (Fig. lA, compare lanes 2 and 5 with 3 and 4). The 
experiment was repeated to confirm that BGLF4, and not BGLF5, 
reproducibly up-regulated expression of FR3 when expressed 
alone or together with BGLF5 (Fig. SI). 

Since FR3 is a late protein, we examined the possibility that lack 
of BGLF4 might disrupt viral DNA replication and thus 
compromise expression of FR3. To assess the effect of BGLF4 
on the extent of viral genome amplification we prepared genomic 
DNA from aliquots of the same cells that were examined for 
protein expression. The relative concentration of viral DNA was 
assessed using primers specific to the upstream region of oriLyt 
and quantitative polymerase chain reaction (qPCR). We found 
that expression of BGLF4 slightly increased viral DNA replication 
by 1 .8-fold relative to cells transfected with ZEBRA alone or with 
ZEBRA plus BGLF4 and BGLF5 (Fig. IB, lanes compare lane 3 
with 2 and 4). Over-expression of ZEBRA and BGLF5 reduced 
the extent of EBV genome amplification by 50% relative to 
ZEBRA alone (compare lanes 2 and 5). Thus lack of BGLF4 lead 
to a modest reduction (l.B-fold) in viral DNA replication but 
caused a more pronounced defect (8.2-fold) on synthesis of FR3 
late protein. 

The kinase activity of BGLF4 is necessary for activation of 
late gene expression but not viral DNA replication 

Lysine 102 in BGLF4 corresponds to a conserved lysine present 
in the catalytic domain of protein kinases [53]. Substitution of 
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Figure 1. BGLF4 markedly enhances expression of the BFRF3 late protein. The experiment Illustrates the effect of Introducing BGLF4 and 
BGLF5 on expression of EBV lytic proteins and viral DNA replication In delta G4/G5 cells. (A) Western blot analysis of EBV early and late proteins 
expressed 48 h after transfectlon of the Indicated plasmlds. (B) Quantitative PGR measuring amplification of the EBV genome from the same samples 
as In panel A. Z, ZEBRA; FR3, BFRF3 minor capsid protein (late); BMRF1, DNA polymerase processlvlty factor (early). 
dol:1 0.1 371 /journal.ppat.1 004307.g001 
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lysine 102 to isoleucine abolishes the kinase activity of BGLF4 
[20]. To determine whether the kinase activity of BGLF4 is 
required for its effects on viral DNA replication and synthesis of 
FR3, we expressed BGLF4 or BGLF4(K102I) together with 
ZEBRA in delta G4/G5 cells and harvested the cells after 48 h. 
We did not detect any difference in viral DNA repKcation as a 
result of expressing the kinase-active or -inactive forms of BGLF4. 
Both wild-type and mutant BGLF4 equally enhanced viral DNA 
replication by 2.6-fold relative to ZEBRA alone (Fig. S2B, lanes 2, 
4, and 6). However a clear difference was observed between 
BGLF4 and BGLF4(K102I) in promoting expression of the FR3 
late gene (Fig. S2A). Wild-type BGLF4 increased expression of 
FR3 by 9-fold whereas the kinase dead mutant of BGLF4 
enhanced the level of FR3 by 1.66-fold. Wild-type and mutant 
forms of BGLF4 protein were equally expressed. The level of 
BMRF 1 protein did not change when the kinase activity of BGLF4 
was disrupted (Fig. S2A), but the hyper-phosphorylated form of 
BMRFl was only detected in cells expressing active BGLF4 and 
not the kinase dead mutant. These findings show that the kinase 
activity of BGLF4 is necessary for enhanced expression of FR3 but 
is not required for the low level of stimulation of viral DNA 
replication by BGLF4. 

Increasing concentrations of wild-type BGLF4, but not 
the kinase dead mutant, up-regulate expression of the 
late FR3 protein independent of effects on viral DNA 
replication 

The previous experiment suggested that the kinase activity of 
BGLF4 is necessary for up-regulating expression of the late FR3 
protein. To further examine this point, we expressed increasing 
concentrations of wild-type or mutant BGLF4 in the presence of a 
constant amount of ZEBRA. After 48 h, transfected delta G4/G5 
cells were harvested and the levels of FR3, BMRFl, BGLF4 and 
ZEBRA proteins were assessed using Western blot analysis. By 
comparing the effect of the two forms of BGLF4, we found that as 
we increased the level of wild-type BGLF4 there was a progressive 
increase in the amount of FR3 protein expressed. The increase in 
the amount of BGLF4 protein expressed was associated with an 
increase in the level of hyper-phosphorylated form of BMRFl 
(Fig. 2A). On the contrary, in cells expressing increasing levels of 
the kinase inactive form of BGLF4, expression of FR3 r(;mained 
constant at a low level and BMRFl was not hyperphosphorylated. 

Using qPCR we studied the effect of expressing various 
concentrations of BGLF4 and BGLF4(K102I) on viral DNA 
replication; we prepared DNA from the same cells used in 
figure 2 A. We only detected a change in the amount of viral DNA 
when we transfected 6 |.lg of wild-type BGLF4. At this input, 
replication increased by 1 .6-fold relative to expression of ZEBRA 
alone (compare lanes 2 and 5). Findings from Figures S2 and 2 
strongly indicate that the effect of BGLF4 on late gene expression 
is kinase dependent and is separate from its effect on viral DNA 
replication. 

BGLF4 regulates late genes transcripts 

In the next experiment we asked whether the effect of BGLF4 
on expression of FR3 could be detected at the transcript level, and 
whether BGLF4 regulated mRNA expression of other late genes. 
We prepared total RNA from fractions of the same samples that 
were previously used to generate figure S2 and employed 
quantitative reverse transcriptase polymerase chain reaction 
(RT-qPCR) to assess the level of five lytic transcripts. Four of 
these EBV transcripts, encoding the major and minor capsid 
proteins (BcLFl and BFRF3), the major glycoprotein gp350/220 



(BLLFl), and a tegument protein (BLRF2), have late kinetics and 
one, encoding BMRFl, has early kinetics. When co-expressed with 
ZEBRA, wHd type BGLF4 or BGLF4(K102I) enhanced the level 
of BMRFl by 2.17 and 1.8-fold, respectively (Fig. 3). This result 
suggested that BGLF4 had a limited stimulatory effect on 
expression of the early BMRF 1 transcript in a manner indepen- 
dent of its kinase activity. However, the kinase activity of BGLF4 
was necessary for efficient expression of the four late transcripts 
examined. The level of all four late transcripts was less in delta 
G4/G5 cells expressing ZEBRA alone or ZEBRA plus G4(K102I) 
compared to co-expression of ZEBRA and wild-type BGLF4 
(Fig. 3). Expression of late transcripts was reduced by an average 
of 7.9-fold (range 9.7 to 5.4) in cells expressing ZEBRA alone and 
a mean of 8.4-fold (range 11.6 to 5.8) in cells expressing ZEBRA 
plus the kinase dead mutant. These results indicate that the kinase 
activity of BGLF4 markedly up-regulates expression of late but not 
early transcripts. 

Knockdown of BGLF4 abolished expression of the late 
BFRF3 protein 

To eliminate the possibility that the markedly enhanced 
synthesis of EBV late products in delta G4/G5 cells following 
ectopic expression of BGLF4 was due to its overexpression, we 
used siRNA to knockdown endogenous BGLF4. For these 
experiments we used BZKO cells, which are 293 cells containing 
an EBV bacmid that lacks a functional gene for ZEBRA, to study 
the effect of silencing BGLF4 on viral DNA replication and late 
gene expression. Transfection of ZEBRA (2 |J.g) induced the lytic 
cycle and activated expression of endogenous BGLF4 (Fig. 4A, 
lane 2). To knock down expression of BGLF4 we co-transfected 
four concentrations, 2, 5, 10 and 15 nM, of siRNA specific to the 
viral protein kinase (siG4) (Fig. 4A, lanes 3 to 6). Cells were 
harvested after 48 h and expression of the FR3 protein was 
assessed using Western blot analysis. We found that all four 
concentrations of siG4 abolished (-xpression of BGLF4. The 
reduction in the level of the late FR3 protein correlated with the 
amount of siG4 used to knockdown expression of BGLF4. This 
reduction in the level of the FR3 protein ranged between 9.8- to 
12.3-fold relative to lytic cells without siG4 (Fig. 4A, lanes 3 to 6). 
We used qPCR to assess the effect of abolishing (expression of 
BGLF4 on replication of the EBV genome. Knocking down of 
endogenous BGLF4 reduced EBV DNA levels by 2.1- to 2.4-fold 
(Fig. 4B, compare lane 2 with lanes 3 to 6). The effect of siRNA to 
BGLF4 on replication was similar in magnitude to that observed 
in previous experiments using delta G4/G5 cells in which 
expression of either wUd-t^pe or the kinase dead mutant of 
BGLF4 enhanced expression of BMRFl and viral DNA replica- 
tion by ~2-fold (Figs. 1, S2 and 2). 

To attempt to segregate the effect of BGLF4 on viral DNA 
replication from that on late gene expression we gradually 
increased the level of BGLF4 expressed by titrating the amount 
of siG4 (15, 10, 5 and 2 nM) co-transfected together with fixed 
concentrations of plasmids encoding ZEBRA (2 ng) and BGLF4 
(3 ng). At the highest concentration of siG4 (15 nM), the level of 
the FR3 protein was markedly reduced by 7.5-fold compared to 
cells without siG4 (Fig. 4A, compare lanes 7 and 8). However, the 
level of viral DNA in the same cells was not compromised as a 
result of knocking down BGLF4; on the contrary- the amount of 
viral DNA present was slightly higher relative to cells expressing 
ZEBRA alone or ZEBRA plus BGLF4 (Fig. 4B, lanes 8-11). As 
the level of BGLF4 expressed increased due to decreasing the 
amount of siG4 transfected, there was a proportional increase in 
the amount of the FR3 late protein, but the level of viral DNA 
replication was unaffected. These results establish a direct 
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Figure 2. Increasing concentrations of wild-type BGLF4, but not the kinase dead mutant, enhances late gene expression. Delta G4/G5 
cells were transfected with ZEBRA alone or together with 1, 3 and 6 |j.g of wild-type or kinase dead mutant BGLF4. Cells harvested after 48 h were 
divided into two aliquots to study expression of lytic proteins and viral replication. (A) Western blot analysis to determine the expression levels of 
BMRFl, BGLF4, ZEBRA and FR3 proteins using specific antibodies. (B) Quantitative measurement of the extent of viral DNA replication using real time 
PCR. 

doi:1 0.1 371 /journal.ppat.1 004307.g002 
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Figure 3. Expression of BGLF4 significantly up-regulates tKie level of late transcripts. The abundance of five EBV lytic transcripts was 
assessed, using quantitative RT-PCR, in delta G4/G5 cells transfected with the following plasmids: empty vector (CIVIV), ZEBRA (Z), BGLF4, ZEBRA plus 
BGLF4, mutant BGLF4(K102I), and ZEBRA plus mutant BGLF4(K102I). BIVIRFI encodes an early transcript encoding the DNA polymerase processivity 
factor. BFRF3, BcLFl, BLLF1 and BLRF2 are late transcripts that encode the minor viral capsid, the major capsid protein, gp350/220 major 
glycoprotein, and a tegument protein, respectively. The level of each of the five transcripts was normalized to empty vector control. 
doi:1 0.1 371/journal.ppat.1 004307.g003 



relationship between the level of BGLF4 and expression of the late 
FR3 protein independent of any effect of BGLF4 on viral DNA 
replication. 

Silent mutations in BGLF4 suppress the effect of siG4 on 
late gene expression 

Specificity of siRNA represents a major concern in this widely 
used approach to down-regulate gene expression. To determine 
whether the effects of siG4 on late gene expression were a result of 
down-regulating BGLF4 or an oflF-target effect, we generated silent 
mutations in the hglf4 gene that abolished complementarity 
between the BGLF4 transcript and siG4 without altering the 
amino acid sequence. Transfection of siG4 is expected to eliminate 
expression of endogenous BGLF4 and lat(^ genes. However, 
providing the siG4-resistant BGLF4 ectopically should restore 
expression of late genes, only if the effect of siG4 is target specific. 
Two forms of the siG4-resistant BGLF4 were generated, kinase 
active and kinase dead, referred to as rG4 and rG4(K102I), 
respectively. To study the effect of siG4 on expression of rG4 and 
rG4(K102I), we transfected 2089 cells, 293 cells containing wild- 
type EBV bacmid, with two concentrations of siG4 (10 and 
20 nM) together with expression vectors of ZEBRA and rG4 or 
rG4(K102I). Expression of ZEBRA by itself in 2089 cells triggered 
the lytic program including expression of the endogenous BGLF4 
(Fig. 5A, lane 2). Both concentrations of siG4 were sufficient to 
abolish expression of ectopic and endogenous BGLF4 in cells co- 
transfected with ZEBRA and BGLF4 (Fig. 5A, lanes 4 and 5). 
Expression of rG4 and rG4(K102I) was not affected by the same 
concentrations of siG4 (Fig. 5 A, lanes 7, 8, 10 and 11). This 
outcome allowed us to assess the specificity of siG4 by restoring 
expression of BGLF4 and studying its effect on late gene 
expression. If siG4 non-specificaUy abolished expression of other 
proteins that are crucial for late gene expression, then restoring 
expression of BGLF4 by rG4 would not be sufficient to rescue 
expression of FR3. We found that expression of rG4 and ZEBRA 
in the presence of siG4 restored expression of FR3 to normal levels 
(lanes 7 and 8). However, providing rG4(Kl()2I) instead of rG4, 
failed to restore the level of FR3 protein (lanes 10 and 11). These 
findings showed that down-regulation of FR3 in cells transfected 
with siG4 is due to abolishing expression of BGLF4 with intact 
kinase activity. Lane 9 shows that cells expressing the transfected 
kinase dead protein and endogenous BGLF4 are permissive for 
expression of FR3. While rG4 restored late gene expression, 
neither rG4 nor rG4(K102I) had any significant effect on viral 
DNA rephcation (Fig. 5B). 

We studied the effect of silencing BGLF4 on the level of eight 
lytic transcripts in 2089 cells expressing ZEBRA either alone or 
together with rG4 or rG4(K102I) (Fig. 6). Two of these transcripts, 
BRLFl and BMRFl, are synthesized prior to viral DNA 
replication. The hrlfl gene encodes the Rta protein, a transcrip- 
tion activator and a replication protein, and the hmrfl gene codes 
for the polymerase processivity factor. Both proteins play vital 
roles in the process of viral DNA replication [54,55]. The other six 
transcripts, with late kinetics, encode the following proteins: major 
glycoprotein gp350/220 (BLLFl), major capsid protein (BcLFl), 
tegument protein (BLRF2), triplex capsid protein (BDLFl), 
scaffold protein (BdRFl), and minor capsid protein (BFRF3). 



Knockdown of endogenous and transfected BGLF4 by siG4 
reduced the level of late transcripts by a median fold-change of 
9.1. Conversely, expression of rG4 restored the amount of late 
transcripts to a level equivalent to that observed in cells transfected 
with ZEBRA alone. However, expression of rG4(K102I) failed to 
up-regulate the level of all six late transcripts and resulted in a 
median reduction of 8.6-fold relative to ZEBRA alone. These 
results demonstrate a specific role for the kinase activity of BGLF4 
in regulating the level of several late viral transcripts. 

Analysis of the EBV transcriptome during lytic induction 
reveals the general importance of BGLF4 in regulation of 
late gene expression 

To study the effect of BGLF4 on the EBV transcriptome during 
lytic induction, we performed deep sequencing analysis on RNA 
purified from 2089 cells transfected with empty vector, ZEBRA or 
ZEBRA plus siG4. To demonstrate that siG4 had silenced BGLF4 
in this particular experiment we assessed the levels of BGLF4 and 
FR3 proteins by Western blot. As expected, ZEBRA induced 
expression of BGLF4 and FR3; addition of siG4 abolished 
expression of BGLF4 and the late FR3 protein (data not shown). 
We used the Expectation-Maximization (EM) algorithm in RSEM 
to estimate gene expression levels and EBSeq to identify 
differentially expressed genes [56,57]. Figure 7 compares the 
number of normalized reads corresponding to each viral transcript 
expressed in the absence and presence of BGLF4. The change in 
expression is represented as a log base 2-fold. Red bars designate 
late genes; blue bars designate latent and early kinetic classes of 
lytic genes. Determination of the kinetic classes of EBV genes was 
based on viral gene classification reported by Yuan et al [58] . Viral 
transcripts with a negative fold-change represent EBV genes 
whose expression is BGLF4-dependent, and vice versa. The most 
striking result was that silencing of BGLF4 down-regulated the 
level of thirty-one late transcripts. The level of only three late 
transcripts, encoding BVRFl, BILFl and BMRF2, increased by 
1.2-, 1.8-, and 3.2-fold, respectively (Fig. 7). BVRFl is a tegument 
protein and is a homolog of HSV UL25 (cork), BILFl is a viraUy 
encoded G protein coupled receptor, and BMRF2 is a 
transmembrane envelope protein [59,60,61]. siRNA to BGLF4 
did not reduce the level of transcripts encoding comijoiicnts of the 
replication machinery. Silencing of BGLF4 impaired expression of 
two lytic transcripts of unknown kinetic class, namely, BGLF3.5 
and BGLF3 (Fig. 7) [58]. The reduction in the level of these two 
transcripts in lytic cells expressing BGLF4 compared to cells 
without BGLF4 was 11.6-fold for BGLF3.5 and 7.3-fold for 
BGLF3. The functions of BGLF3.5 and BGLF3 in the EBV lytic 
cycle remain to be characterized. However, orf34, a homologue of 
BGLF3 encoded by murine gamma herpesvirus-68 (MHV-68), is 
an early gene that was found to be essential for late gene 
expression and not viral DNA replication [4]. These findings 
demonstrated that BGLF4 globally enhanced expression of EBV 
late genes. 

Knockdown of BGLF4 down-regulates expression of 
BGLF3 and BGLF3.5 

RNA-seq analysis showed that silencing of BGLF4 reduced 
expression of BGLF3 and BGLF3.5. We sought to validate the 
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Figure 4. siRNA specific to BGLF4 abolishes expression of tKie late FR3 protein. (A) Immunoblot of BZKO cell lysates expressing ZEBRA 
alone (lanes 2 to 6) or ZEBRA and BGLF4 (lanes 7 to 1 1 ) together with increasing concentrations (2 to 1 5 nIVI) of siRNA complimentary to BGLF4 (siG4). 
The levels of BGLF4 and late BFRF3 (FR3) proteins were determined by densitometry. (B) Fold change in viral lytic DNA replication as a result of 
transfecting siRNA specific to BGLF4 (siG4). Total DNA was prepared from the same cells used in the immunoblot in panel A. 
dor:l 0.1 371/journal.ppat.1 004307.g004 



effect of siG4 on the level of transcripts encoding BGLF3 and 
BGLF3.5 using RT-qPCR. The data presented in figure S3 
represent the mean of two biological replicates. Similar to our 
observations with the RNA-seq analysis, we found that expression 
of the bglf3 and hglf3.5 genes was reduced in cells lacking BGLF4 
by 6-fold and 22-fold, respectively. 

The effect of siG4 on BGLF3 and BGLF3.5 could be attributed 
to two different mechanisms: siRNA to BGLF4 simultaneously 
reduces the BGLF3 and BGLF3.5 transcripts; BGLF4 regulates 
synthesis or stability of BGLF3 and BGLF3.5 mRNAs. The hglf4, 
bglf3.5 and bglfS genes are located in a transcription unit with a 
single poly(A) signal (Fig. 8D). Previous reports suggested that 
transcripts synthesized from this region are co-terminal [33,62]. 
To assess the number and size of transcripts expressed from this 
locus we purified RNA from 2089 cells transfectcd with ZEBRA. 
Using Northern blot analysis and probes directed to the unique 
sequences in BGLF4 and BGLF3.5 we detected five RNAs that 
were 4.3, 3.5, 3.3, 2.2 and 1.8 kb long (Fig. 8A and C). Based on 
the lengths of the genes and the distribution of TATA boxes 
present in this locus, the 4.3 kb RNA correlates with transcript A 
(BGLF3-BGLF3.5-BGLF4-BGLF5-BBLF1). siRNA to BGLF3 
abolished the 4.3 kb transcript without affecting the level of the 
other transcripts (Fig. 8C). This result demonstrates that the 
4.3 kb RNA is the only transcript that contains the BGLF3 
sequence, which is consistent with the composition of transcript 
A. The 3.5 kb RNA corresponds with transcript B (BGLF3.5- 
BGLF4-BGLF5-BBLF1); the 3.2 kb RNA with transcript C 
(BGLF4-BGLF5-BBLF1), and the 1.8 kb RNA with transcript D 
(BGLF5-BBLF1) (Fig. 8D). The 1.8 kb RNA was detected by 
the BGLF4 probe presumably as a result of containing BGLF4 
sequence in the 5' UTR of RNA species D (Fig. 8A). A 2.2 kb 
RNA was only detected using the BGLF4 probe (Fig. 8A); its 
size does not correlate with any of the predicted RNA 
transcripts. 

siRNA to BGLF4 aboUshed four RNA species: 4.3, 3.5, 3.3 
and 2.2 kb. This outcome conforms to the predicted composi- 
tion of transcripts A, B, and C, in which all three transcripts 
contained the BGLF4 sequence (Fig. 8D). Thus, targeting 
BGLF4 with siG4 silences expression of BGLF3.5 and BGLF3 
as well. Transfection of siG4 resulted in two new RNA species 
that were 1.6 and 3 kb long and were only detected by the 
BGLF4 probe (Fig. 8 A). The appearance of these two RNAs is 
consistent with the capacity of siG4 to trigger cleavage of the 
BGLF4 containing transcripts. 

To determine whether BGLF4 has the capacity to increase the 
level of BGLF3 mRNA either by up-regulating its expression or 
enhancing the stability of its mRNA, we compared the level of 
BGLF3 mRNA in 2089 cells transfected vnth ZEBRA, ZEBRA 
plus siG4 or ZEBRA plus siG4 and rG4 (siG4 resistant BGLF4). 
Addition of siG4 reduced the level of BGLF4 and BGLF3 
transcripts expressed from the endogenous viral genome. Ectopic 
expression of rG4 increased the level of BGLF3 transcript by 4- 
fold relative to cells transfected with ZEBR,^ plus siG4 (Fig. S4). 

These results suggest that siG4 reduces the level of BGLF3 
mRNA by two different mechanisms; it directly silences 
expression of BGLF3, and it knocks down expression of BGLF4, 
a protein that augments expression or stability of the BGLF3 
transcript. 



Knockdown of BGLF3 abolished expression of EBV late 
genes 

To investigate the possibility that BGLF4 might exert its effects 
on late gene expression by controlling expression of a second tier 
of late gene regulators we studied the effect of siRNA to BGLF3 
and BGLF3.5 on synthesis of late products. 2089 cells transfected 
with two different concentrations of ZEBRA expression vector (1 
and 2 |J.g DNA) in the presence of siRNA were harvested after 
48 h. Cell lysates were analyzed using Western blot assays. The 
experiment illustrated in figure 9A is a representative of two 
biological replicates. Knockdown of BGLF3 significantly reduced 
expression of the late FR3 protein by an average of 9.5-fold, but 
siRNA to BGLF3 (siG3) had no effect on expression of ZEBRA, 
BGLF4, or EA-D (Fig. 8B and 9A). Silencing of BGLF3.5 reduced 
FR3 expression by an average of 2.7-fold. However, siRNA to 
BGLF3.5 (siG3.5) slightly reduced expression of BGLF4 but not 
ZEBRA or EA-D. To assess the effect of both siRNAs on viral 
DNA replication we purified DNA from the same cells and 
examined the level of viral DNA using qPCR. Neither of the two 
siRNAs reduced the level of viral DNA in cells expressing ZEBRA 
(Fig. 9B). We used RT-qPCR to demonstrate that the siRNA to 
BGLF3 reduced the level of the endogenous BGLF3 transcript by 
2.7-fold, but had no effect on the level of BMRFl mRNA (Fig. S5). 
Expression of four other late transcripts, BDLFl, BdRFl, BLLFl, 
and BcLFl, was reduced as a result of silencing BGLF3 (Fig. S6). 
These results demonstrate that BGLF3 is necessary for progress of 
the EBV lytic cycle into the late phase. 

siG3-resistant BGLF3 (rG3) rescued expression of late 
genes 

To demonstrate that the effect of siG3 on late gene expression 
was specific to knockdown of the bglf3 gene, we inserted silent 
mutations in the region of the BGLF3 mRNA that is recognized 
by siG3. These mutations disrupt the complementarity between 
siG3 and the BGLF3 transcript without altering the amino acid 
sequence of the protein. Insertion of these mutations resulted in a 
form of BGLF3, referred to as rG3, which is resistant to siG3. Co- 
transfection of ZEBRA plus siG3 in 2089 cells reduced the level of 
FR3 by 7-fold relative to cells solely transfected with ZEBRA 
(Fig. lOA, lanes 2 and 3). Ectopic expression of rG3 in cells 
transfected with ZEBRA and siG3 restored expression of the late 
FR3 protein to a level equivalent to that observed in cells 
expressing ZEBRA alone (Fig. lOA, compare lanes 2 and 4). The 
experiment was repeated twice and similar results were obtained. 
To assess the effect of rG3 on expression of other late genes, we 
purified total RNA from the same set of samples and examined 
expression of BMRFl, an early gene, and two late genes: BcLFl, 
and BdRFl. Using RT-qPCR, we found that siG3 reduced die 
level of the two late transcripts but had no effect on the level of the 
early BMRFl transcript; co-transfection of rG3 rescued expression 
of BcLFl and BdRFl (Fig. lOB). These results demonstrate that 
BGLF3 is indispensable for expression of late genes. 

Both BGLF4 and BGLF3 are indispensable for late gene 
expression 

To determine whether BGLF3 or BGLF3.5 can substitute for 
the role of BGLF4 in activation of late gene expression, we 
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Figure 5. Silent mutations in BGLF4 override tKie effect of siG4 on late gene expression. (A) Western blot analysis of 2089 cell lysates 
transfected with the indicated expression vectors. Two concentrations of siRNA specific to BGLF4 (siG4), 10 and 20 nM, were co-transfected with 
expression vectors encoding ZEBRA and wild-type BGLF4 or kinase dead BGLF4(K1 021) (lanes 4, 5, 10 and 1 1 ). To demonstrate the specificity of siG4, 
silent mutations were introduced into the kinase active (rG4) and kinase dead (rG4(K1 021)) bglf4 genes that rendered the BGLF4 transcript resistant to 
siG4. Lysates of cells harvested 48 h after transfection were analyzed by immunoblotting using antibodies specific to BGLF4, ZEBRA and the late 
BFRF3 protein (FR3). (B) Quantitative PCR to measure the extent of EBV genome amplification in each condition. 
dor:1 0.1 371 /journal.ppat.1 004307.g005 



knocked down BGLF4 by transfecting siG4 and expressed BGLF3 
and BGLF3.5, separately and together, using expression vectors. 
We found that provision of BGLF3 and BGLF3.5 was not 
sufficient to trigger expression of the late FR3 protein in absence of 
BGLF4 (Fig. S7). Similarly, ectopic expression of BGLF3 and 
BGLF3.5 did not restore late gene expression in cells treated with 
phosphonoacetic acid (PAA), an inhibitor of viral DNA replication 
(Fig. S7, lanes 8 and 9). These restdts demonstrate that BGLF4, 
BGLF3 and viral DNA replication are three necessary compo- 
nents for progress of the EBV lytic cycle into the late phase. 

Discussion 

Classification of lytic viral gene expression into pre- and post- 
replication temporal phases is applicable to the life cycles of all 
herpesviruses. While kinetics of expression of pre-replication genes 
has been extensively studied, little is known about the order of 
events that take place prior to expression of late genes. Our 
findings reveal a previously unknown transitional step between EB 
viral DNA replication and late gene expression. We demonstrate 
that the kinase function of BGLF4 is necessary for optimal 
expression of late genes. Inactivating the kinase activity of BGLF4 
markedly impairs late gene expression without impeding viral 
DNA replication; the kinase dead mutant supported viral DNA 
replication to the same extent as wild-type BGLF4 (Fig. S2, 2, and 
5). Disrupting the kinase activity or the expression of BGLF4 
selectively reduced the level of late transcripts (Fig. 3, 6 and 7). 
Analysis of the EBV transcriptome revealed that expression of 31 
late genes was reduced in the absence of BGLF4 (Fig. 7). Silencing 
of BGLF4 did not reduce the level of most early transcripts 
including those encoding replication proteins. Using siRNA to 
compare the expression pattern of EBV genes with and without 
BGLF4 we identified two non-late genes, bglf3 and hglf3.5, whose 
expression was significantly reduced as a result of silencing BGLF4 
(Fig. 7, 8 and S3). By knockdown experiments, we found that 
BGLF3 is an independent regulator of late g('ne expression. siRNA 
to BGLF3 markedly reduced the level of FR3, a canonical late 
protein, and the transcripts of several other late genes without 
affecting endogenous expression of BGLF4 or viral DNA 
replication (Fig. 8B, 9A, 10 and 86). Expression of siRNA resistant 
forms of BGLF3 and BGLF4 provided evidence that each protein 
plays an independent and indispensable role in regulation of late 
gene expression (Fig. 5, 6, and 10). In summary, our results 
identify a control locus composed of BGLF4 and BGLF3 that 
regulates synthesis of EBV encoded structural proteins (Fig. 11). 

Role of BGLF4 in regulating viral DNA replication and late 
gene expression 

To establish a distinct regulatory role of BGLF4 in activation of 
late gene expression it is imperative to exclude that this role results 
from an effect on viral DNA replication. Several groups have 
assessed the effect of knocking down or knocking out BGLF4 on 
the process of EBV genome amplification. Two separate research 
groups reported that knocking out the hglf4 gene did not affect the 
level of viral DNA synthesized during the lytic cycle [28,32]; three 
other groups found that lack of BGLF4 reduced viral DNA 



rephcation by 1.4- to 2-fold [31,44,47]. We found that ectopic 
expression of BGLF4 in BGLF4 nuU cells enhanced viral DNA 
replication by 1.6 to 2.6-fold (Fig. 1, S2, and 2) in agreement with 
the second group of studies. These slight differences in the effect of 
BGLF4 on viral DNA replication could be attributed to the cell 
background or the approach used to knockout the hglf4 gene. In 
our knockdown experiments we found that silencing of BGLF4 in 
BZKO ceUs reduced viral DNA replication by 2-fold and in 2089 
ceUs by 1.1 -fold (Fig. 4 and 5). The effect of BGLF4 on viral DNA 
replication was independent of its kinase activity. In delta G4/ G5 
cells, expression of either wild-type BGLF4 or the kinase inactive 
mutant enhanced rephcation to same extent (Fig. 2 and 82). Thus, 
we conclude that the kinase activity of BGLF4 is required for its 
effect on late gene expression but not for its modest enhancement 
of viral DNA rephcation. 

Since expression of late genes is stipulated by the onset of viral 
genome amplification we sought to segregate the effect of BGLF4 
on replication from its effect on expression of late genes. We 
achieved this goal using two difierent approaches: first, we varied 
the concentration of siRNA to BGLF4. We demonstrated that at a 
very low concentration of BGLF4 protein there was significant 
reduction in late gene expression without any significant compro- 
mise in viral DNA replication (Fig. 4, compare lane 7 with lane 8). 
As we increased the concentration of BGLF4, viral DNA 
replication remained unchanged (Fig. 4B) while late gene expres- 
sion increased proportionally (Fig. 4A). Second, we supplied a 
kinase dead form of BGLF4 to cells with a knockout of the hglf4 
gene. Viral DNA replication was rescued to normal levels but 
there remained a marked defect in late gene expression (Fig. 82 
and 2). 

The possibility that BGLF4 possesses additional functions that 
are independent of its kinase activity is supported by the capacity 
of ORF36, a BGLF4 homolog encoded by MHV68, to repress the 
function of class 1 and 2 histone deacetylases (HDACs) [43]. 
Repression of HDAC 1 and 2 by ORF36 is mediated by protein- 
protein interactions, rather than by phosphorylation. Expression of 
either wild-type or kinase inactive ORF36 enhances MHV68 early 
gene expression and viral DNA replication [43]. The effect of 
BGLF4 or BGLF4(K102I) on EB viral DNA replication may also 
be due to the capacity to interact with and disrupt the function of 
HDACl and 2 during the early phase of the lytic cycle. BGLF4 
interacts with HDACl and 2 [43]; the biological significance of 
such interaction for temporal control of the EBV fife cycle has not 
been studied. 

Expression of BGLF4, BGLF3.5 and BGLF3 from the late 
control locus in the EBV genome 

Expression profiling of EBV genes during lytic infection affirms 
that most late genes are subject to regulation by BGLF4 (Fig. 7). 
By comparing the EBV transcriptome during lytic infection in the 
absence and presence of BGLF4 we found that BGLF3 and 
BGLF3.5 represent the two most down-regulated non-late genes in 
cells transfected with siG4 (Fig. 7). The effect of siG4 on the level 
of BGLF3 and BGLF3.5 mRNAs could be attributed to two 
different mechanisms. First, in addition to targeting the BGLF4 
mRNA, siG4 might concurrendy knock down the transcripts 
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Figure 6. Expression of siG4-resistant BGLF4 restores synthesis of EBV late transcripts. The effect of siG4 on the level of two early and six 
late viral transcripts was analyzed by RT-qPCR. The early transcripts are BIV1RF1 (polymerase processivity factor) and BRLF1 (transcription factor). The 
late transcripts are BLLF1 (gp350/gp220 major glycoprotein), BFRF3 (minor capsid protein), BLRF2 (tegument protein), BDLF1 (triplex capsid protein), 
BdRFI (scaffold protein), and BcLFI (major capsid protein). The relative abundance of each transcript was calculated using the AAC7 metiiod. 
doi:l 0.1 371/Journal.ppat.1 004307.g006 



encoding BGLF3 and BGLF3.5. hglfS, hglf3.5 and bglf4 are 
nested within a transcription unit containing five overlapping open 
reading frames that also includes hglfS and hblfl, respectively 



(Fig. 1 1). Due to the presence of a single canonical poly(A) signal 
at the 3 ' end of hblfl , transcription from this locus is likely to result 
in co-terminal transcripts (Fig. 8D) [33,62], unless cryptic variants 
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Figure 7. Impact of knockdown of BGLF4 on the EBV transcriptome. RNA-seq analysis compared EBV transcripts expressed in 2089 cells 
transfected with ZEBRA versus ZEBRA plus siG4. RNA was prepared at 48 h after transfection. The plotted graph was generated by dividing the 
number of normalized reads of a particular transcript in the presence of siG4 by the number of reads obtained in the absence of siG4. This ratio 
represents the effect of BGLF4 on the abundance of a particular transcript and is presented as a log base 2-fold. A negative fold-change indicated 
dependency on BGLF4 for expression. Late genes illustrated in red bars while latent and other lytic genes were illustrated in blue bars. 
doi:10.1371/journal.ppat.1004307.g007 



of the poly(A) signal exist at the 3 'end of some genes [63]. In 
figure 8A, knockdown of BGLF4 abolished all the RNA species 
that were predicted to contain BGLF4 sequence. This result is in 
favor of co-terminal expression of transcripts encoded by the late 
control locus. Second, BGLF4 might regulate the level of the 
BGLF3 and BGLF3.5 mRNAs. In Fig. S4, siG4 down-regulated 
the level of BGLF3 mRNA; provision of rG4 overcame the effect 
of siG4 and increased the level of the endogenous BGLF3 
transcript by 4-fold. BGLF4 might either enhance the activity of 
the BGLF3 promoter or enhance the stability of the BGLF3 
mRNA. 

Silencing of BGLF3 had no discernable effect on the protein 
level of BGLF4 (Fig. 8B and 9A). A possible explanation of this 
result is that BGLF3 is expressed at a low level during the late 
phase of the lytic cycle. Northern blot analysis revealed that the 
4.3 kb transcript, the only RNA containing the BGLF3 sequence, 
was considerably less abundant relative to the other RNA species 
synthesized from this locus (Fig. 8A and C). Therefore, knock- 
down of the 4.3 kb RNA (Fig. 8D, transcript A) presumably had 
more impact on the level of BGLF3 relative to the overall level of 
BGLF4 containing mRNAs (Fig. 8C and D). In two independent 
experiments, abolishing the BGLF3 containing transcript by siG3 
reduced the abundance of the FR3 protein without affecting the 
level of the BGLF4 protein (Fig. 8B and 9A). 

Contrary to siG3, siRNA to BGLF3.5 (siG3.5) reduced the level 
of the BGLF4 protein by 3-fold (Fig. 9A, lanes 3 and 7). BGLF4 
has three transcriptional start sites; two sites mapped to regions 
within the hglf3.5 coding sequence and a third site mapped 
upstream of the translation start codon of hglf3.5. Reduction in 
the level of BGLF4 protein by siG3.5 is likely due to the ability of 
siG3.5 to target transcripts that encode both BGLF3.5 and BGLF4 
(Fig. 1 1). In KSHV, ORF36, the homolog of BGLF4, is robusdy 
translated as a downstream cistron from a polycistronic transcript 
that initiates with ORF35, the KSHV homolog of BGLF3.5 
[64,65]. 



Both BGLF4 and BGLF3 are independently necessary for 
late gene expression 

The kinase activity of BGLF4 and not the mere synthesis of the 
BGLF4 transcript or protein is necessary for stimulation of late 
gene expression (Fig. 2, S2, 3, 5, and 6). In knockout and 
knockdown experiments of BGLF4, expression of the kinase 
active, and not the kinase dead, form of BGLF4 activated synthesis 
of late transcripts. Attempts to complement the lack of BGLF4 
with expression of BGLF3, BGLF3.5 or both did not restore 
expression of late genes (Fig. S7). Therefore, the role of BGLF4 in 
late gene expression is not limited to regulating the level of 
BGLF3; BGLF4 has an unknown but indispensable function in the 
mechanism regulating synthesis of late products. BGLF3 is an 
ortholog of the MHV68-encoded ORF34 [4]. Mutant MHV68 
virus that lacks the orf34 gene underwent viral genome 
amplification but failed to express late products. Here, we found 
that the function of MHV-68 ORF34 in inducing late gene 
expression is conserved in EBV BGLF3. Silencing expression of 
BGLF3 using siRNAs abolished synthesis of the FR3 late protein 
and the transcripts of several late genes without affecting the level 
of BGLF4 or viral DNA replication (Fig. 8, 9, and S6). Providing a 
form of BGLF3 that is resistant to siG3 annulled the effect of the 
siRNA on expression of late genes. The exact role of BGLF3 in 
regulation of EBV late gene expression has not been established. 
Knockout of ORF34, the murine homolog of BGLF3, abolished 
recruitment of RNAPII to late promoters. In a related observation, 
we found that BGLF3 interacts with the C-terminal domain of 
RPBl, the large catalytic subunit of RNAPII (data not shown). 
These results represent the first demonstration of the indispensable 
roles of BGLF4 and BGLF3 in regulation of EBV late gene 
expression. However, ectopic expression of BGLF4 or BGLF3 
does not complement the lack of viral DNA replication (Fig. S7, 
lanes 8 and 9, and data not shown). Therefore, viral DNA 
replication as well as the function of BGLF4 and BGLF3 represent 
three essential components for stimulation of late gene expression. 
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Figure 8. Knockdown of BGLF3 has no effect on expression of BGLF4. To examine the effect of silencing BGLF3 expression on the mRNA and 
protein levels of BGLF4, we transfected 2089 cells with ZEBRA alone or together with siG3 or siG4. Cells were harvested after 48 h. RNA and protein 
samples were prepared for Northern (A and C) and Western (B) blot analyses. In cells transfected with ZEBRA, we detected four RNA species: 3.5, 3.3, 
2.2 and 1.8 kb using the BGLF4 probe (A). The 4.3 kb was only detected using a probe unique to BGLF3.5 (C). Except for the 4.3 kb, the same RNA 
species were detected when siG3 was co-transfected. However, the 3.5, 3.3 and 2.2 kb RNA species disappeared in the sample co-transfected with 
siG4. These RNAs were replaced with two new species approximately 1.6 and 3 kb in size. Examining protein expression in the same samples 
demonstrated that siG4 and siG3 markedly reduced expression of the late FR3 protein (Fig. 8B). Nonetheless, similar to Fig. 9, siG3 had no effect on 
the level of the BGLF4 protein. The experiment was repeated twice and similar results were obtained. (D) A schematic diagram of the various 
transcripts expected to be synthesized from the late control locus. CMV, empty vector; Z, ZEBRA; siG3, siRNA to BGLF3; siG4, siRNA to BGLF4. The 
asterisks denote probable cleavage products. 
doi:1 0.1 371 /journal.ppat.1 004307.g008 



Our findings have scientific and translational impUcations. The 
findings provide new insight on regulation of late gene expression, 
one of the main puzzles in virology, and emphasize the importance 
of targeting the kinase activity of BGLF4 for development of new 
antiviral drugs. 

Materials and Methods 

Expression vectors 

The vector encoding the ZEBRA protein was prepared as 
previously described [66]. Constructs expressing wild- type BGLF4 
and the kinase dead BGLF4(K102I) were a kind gift of Dr. Mei- 
Ru Chen [16]. BGLF4 forms that are resistant to the 
corresponding siRNA (siG4) were generated by introducing silent 



point mutations using the following mutagenic primer: 5'- 
GTGACCAACATTGATGACATGACGGAGACATTATACG- 
TCAAATTACCTGAAAACATGACGCGCTGTGATCACCT- 
CCCCATTACC-3' and its complementary strand: 5'-GG- 
TAATGGGG AGGTGATC AC AGC GC GTC ATGTTTTC AG- 
GTAATTTGACGTATAATGTCTCCGTCATGTCATCAAT- 
GTTGGTCAC-3'. 

Cell culture and transfection 

2089, Bam Z knockout (BZKO), and delta BGLF4/BGLF5 
(delta G4/G5) are 293 human embryonic kidney (HEK) cells 
stably transfected with bacmids containing wild-type, BZLFl null, 
and BGLF4/BGLF5 null EBV B9,5.8 genomes, respectively 
[47,67,68]. Mutations in wild type EBV bacmid (2089) were 
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Figure 9. BGLF3 regulates late gene expression. A) Western blot analysis of EBV lytic proteins to study the effect of knocking down expression 
of BGLF3 and BGLF3.5. 2089 cells were transfected with two concentrations of the ZEBRA expression plasmid and siRNAs to BGLF3 (siG3) and BGLF3.5 
(siG3.5). Cells were harvested after 48 hours of transfection. The membrane was blotted with antibodies specific to BGLF4, ZEBRA and the late FR3 
protein. The level of p-actin used as a loading control. B) Total DNA was prepared from the same samples and was analyzed by quantitative PCR for 
the level of the upstream region of oriLyt as a marker for EBV replication. The relative DNA concentration was calculated using the standard curve 
method. The experiment is a representation of two biological replicates. 
doi:1 0.1 371 /journal.ppat.1 004307.g009 
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Figure 10. Silent mutations in BGLF3 suppress the effect of siG3 on late gene expression. (A) Western blot analysis of 2089 cells 
transfected with the Indicated molecules. The membrane was blotted with specific antibodies against three EBV lytic proteins: ZEBRA, BMRF1 (E: 
early), and BFRF3 (L: late). FLAG antibody was used to detect the siG3-resistant BGLF3 (rG3). (B) RT-qPCR to assess the capacity of rG4 to rescue 
expression of two late transcripts, BcLFl and BdRFI. PAA, phosphonoacetic acid, an inhibitor of viral DNA replication. 
doi:1 0.1 371 /journal.ppat.1 004307.g01 0 



generated by homologous recombination in which the bzlfl gene 
and the hglf4 gene were replaced with the kanamycin resistance 
gene [47,67]. Disruption of the hglf4 gene abolished expression of 
BGLF4 and BGLF5. The cells were cultured in Dulbecco's 
modified Eagle medium (DMEM) supplemented with 10% fetal 
bovine serum (FBS), and antibiotics (peniciUin-streptomycin at 50 
units/ml and amphotericin B at 1 |ag/ml). Hygromycin B 
(Calbiochem) 100 |J.g/ml was added to the medium to select for 
293 cells containing wild type or mutant EBV bacmids. 
Transfection of eukaryotic plasmids was performed in 25 cm^ 
flasks using 36 |J.l DMRIE-C transfection reagent (Invitrogen) 
mixed with 2 [ig of ZEBRA expression vector to induce the lytic 
cycle, and 3 |J.g of plasmids encoding wild-type or kinase dead 
BGLF4. 

Silencing expression of BGLF4, BGLF3 and BGLF3.5 was 
attained using RNA interference technique. Transfections were 
performed using Lipofectamine 2000 (Invitrogen) following the 
manufacturer's protocol. Several siRNAs were tested to knock- 
down BGLF4, BGLF3, and BGLF3.5. The most specific siRNA 
was then used in knockdown experiments to study the role of each 



of these proteins during the lytic cycle. AH transfections were 
carried out in OPTI-MEM medium. Cells were incubated at 37°C 
in 5% CO2 incubator and harvested after 48 h of transfection. 

Protein detection 

Harvested cells were re-suspended in sodium dodecyl sulfate 
(SDS) sample buffer at 10* cells/ 10 |J.l. Proteins were separated in 
10% SDS-polyacrylamide gels and transferred to nitrocellulose 
membranes (Bio-Rad). The membrane was blotted with specific 
antibodies to cellular and viral proteins. The BGLF4 antibody was 
raised against amino acids 1 to 220 in rabbit. BGLF4(l-220) was 
expressed in Escherichia coli from a pET-22b vector and purified 
using nickel affinity chromatography. The EA-D (BMRFl) 
monoclonal antibody (R3.1) was kindly provided by G. Pearson 
[69]. Anti-ZEBRA and anti-BFRF3 are polyclonal rabbit 
antibodies and were described previously [54]. FLAG-tagged 
proteins were detected with anti-FLAG mouse monoclonal 
antibody (Sigma), fi-actin was detected using a mouse monoclonal 
antibody (Sigma). The antigen-antibody complex was detected by 
autoradiography using '^'I-protein A. 
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Figure 11 . A model depicting regulation of late gene expression by the G-locus. A schematic diagram illustrating the organization of the 
bglf4 (G4), bglfS.S (G3.5) and bglfS (G3) open reading frames in the EBV genome. Four canonical TATA box elements located in this locus are present 
upstream of the bglf3, bglf3.5, bglfS and bblfl genes [33,47,73]. The asterisks denote three mapped transcription start sites of transcripts encoding 
BGLF4. A poly(A) signal is located at the 3' end of the bblfl gene. Solid lines represent the predicted transcripts synthesized from this locus. Our 
findings postulate that both BGLF4 and BGLF3 function to activate expression of late genes. 
doi:10.1371/journal.ppat.1004307.g011 



Quantitative RT-PCR 

RNA was prepared from cells harvested 48 h after transfection 
using Qia-shredder and RNeasy Plus products from Qjagen. The 
concentration of RNA in each sample was determined by 
measuring the optical density at 260 nm. The level of viral 
transcripts was assessed in 100 ng of total RNA using iScript One- 
Step RT-PCR with SYBR green (Bio-Rad) in a total volume of 
25 |J,1. The level of 18S rRNA was measured to normalize for the 
total amount of RNA. Each sample was analyzed in triplicates and 
the fold change in expression was calculated using the AACf 
formula. The sequences of the primers used are available upon 
request. 

Next generation RNA sequencing and data analysis 

Three strand-specific sequencing libraries were produced from 
total RNA of 2089 cells transfected with empty vector, ZEBRA or 
ZEBRA plus siG4 as previously described [70]. The libraries were 
run on HiSeq 2000, generating approximately 86, 47 and 49 
million reads, respectively. The generated reads were single-end 
and each read was 76 bp long. Adapter sequences, empty reads, 
and low-quality sequences were removed. Also, the first 1 6 and last 
4 nucleotides in each read were trimmed using the FASTX toolkit 
(http://hannonlab.cshl.edu/fastx_toolkit/index.html) to remove 
low quality bases. Trimmed reads were mapped to the human 
reference genome (hgl9) with a known transcriptome index 
(UCSC Known Gene annotation) using Tophat v2.0.8 [7 1]. Those 
reads that did not map to the human genome were later mapped 
to the EBV genome (GenBank accession number NC_007605.1) 
with a known transcriptome annotation [62]. We used the 



Expectation-Maximization (EM) algorithm in RSEM [57] with 
Bowtie 2 [72] to estimate gene expression levels. EBSeq within 
RSEM pipeline was used to identify differentially expressed genes 
[56]. 

Replication assay 

Cells were harvested and re-suspended in 400 [l\ of lysis buffer 
containing 50 mM Tris-HCl [pH 8.1], 1% SDS, and 10 mM 
EDTA. The cell lysate was sonicated twice, 10 pulses each, and 
centrifuged for 10 min in a microfuge at full speed. The 
supernatant was diluted 10-fold in dilution buffer containing 
16.7 mM Tris-HCl [pH 8.1], 0.01% SDS, 1.1% Triton X-100, 
167 mM NaCl, and 1.2 mM EDTA. The samples were subjected 
to Proteinase-K digestion (Roche) followed by phenol-chloroform 
extraction to remove proteins. The DNA was purified using 
Qjagen PCR-purification kit. The total DNA content was 
determined by measuring the absorbance at 260 nm. The extent 
of viral genome amplification was quantitated using the IQ^ Sybr 
Green SuperMix kit (Bio-Rad). The sequences of the forward and 
reverse oriLyt primers are: 5'-TCCTCTTTTTGGGGTCT- 
CTG-3' and 5'-CCCTCCTCCTCTCGTTATCC-3'. The rela- 
tive concentration of DNA was calculated based on a standard 
curve constructed from different concentrations of oriLyt. The 
level of viral DNA was normalized to control sample transfected 
with empty vector (CMV). 

Northern blot. Cell pellets were collected at 48 h after 
transfection with CMV, ZEBRA, and siRNAs. Total RNA was 
isolated using the RNeasy Plus kit (QIAGEN) and the resulting 
RNA was quantified via UV spectrometry. Thirty micrograms of 
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RNA per sample was electrophoresed on a 1% agarose-6% 
formaldehyde gel in 20 mM MOPS {morpholinepropanesulfonic 
acid; pH 7). The gel was transferred to a nylon membrane 
(HYBOND-N+; Amersham Biosciences) and hybridized with ''^P- 
labeled probes. The radioactive probes were complementary to 
unique sequences in BGLF4, BGLF3.5, and RNaseP HI. Probes 
were prepared by random priming. The BGLF4 probe was 
generated from a 620-bp PGR product amplified using the 
following primers: 5'-ATGGATGTGAATATGGGTGG-3' and 
5'-TCTGTGAAATCCACCAGGAT-3'. The BGLF3.5 probe is 
a PGR fragment that was amplified using the following primers: 
5'- TAGAAAAGAGAGGGGCTGTG-3' and 5'-GCGGTGT- 
AGGTGGTTGTAGT-3'. The RNaseP HI probe was prepared 
using a 350-bp NcoI-to-PstI fragment of cDNA from the HI RNA 
component of human RNaseP. 

Supporting Information 

Figure SI BGLF4 induces expression of the late FR3 
protein. In Fig. lA, we showed that expression of BGLF4, but 
not BGLF5, in delta G4/G5 cells up-regulated expression of FR3. 
To confirm this result we repeated the experiment in the same cell 
line and found that BGLF4 is necessary for efficient expression of 
the late FR3 protein. Gomplementation with BGLF5 did not up- 
regulate the level of FR3 protein (lane 4). 
(TIF) 

Figure S2 The kinase activity of BGLF4 is necessary for 
its role in activating late gene expression. (A) Expression of 
viral lytic proteins in delta G4/G5 cells transfected with empty 
vector (GMV), ZEBRA (Z), BGLF4, ZEBRA plus BGLF4, kinase 
dead BGLF4(K102I), and ZEBRA plus BGLF4(K102I). The blot 
was sequentially reacted with antibodies specific to BFRF3, 
ZEBRA, BGLF4, and BMRFl. Since BGLF4 and BMRFl 
cannot be fuUy resolved on 10% SDS-PAGE, figure 2A has two 
panels displaying the BGLF4 protein with and without BMRFl. 
(B) Fold change in the level of EBV DNA was determined by 
quantitative PGR in delta G4/G5 cells transfected with the 
indicated expression vectors relative to empty vector. 
(TIF) 

Figure S3 Knockdown of BGLF4 down-regulates expres- 
sion of BGLF3 and BGLF3.5. Quantitative RT-PCR was used 
to measure the level of the BGLF3 (A) and BGLF3.5 (B) transcripts 
in 2089 cells with and without silencing expression of BGLF4. 
Total cellular RNA was harvested at 48 h after transfection. The 
relative level of each transcript was calculated as a fold-change 
relative to empty vector using the AACr method. 
(TIF) 

Figure S4 siG4-resistant BGLF4 (rG4) enhances expres- 
sion of BGLF3. Knockdown of BGLF4 markedly reduced the 
level of the BGLF3 transcript. Two possible interpretations could 
be envisioned: i) BGLF4 regulates expression of BGLF3 and ii) 
BGLF4 and BGLF3 are encoded by the same transcript. To 
determine if ectopic expression of BGLF4 can up-regulate the level 
of the BGLF3 mRNA transcribed from the endogenous viral 
genome, we knocked down the endogenous BGLF4 transcript 
using siG4 and co-transfected a form of BGLF4 that is resistant to 
the siRNA. Using RT-qPGR we found diat sUencing of BGLF4 
reduced the level of BGLF3 by 2.8-fold relative to cells expressing 
ZEBRA. Go-expression of rG4 enhanced expression of BGLF3 by 
4-fold compared to cells transfected with ZEBRA plus siG4. This 



result indicates that BGLF4 has the capacity to induce expression 

of the endogenous BGLF3 transcript. 

(TIF) 

Figure S5 siG3 specifically reduces the level of the 
BGLF3 transcript. In Fig. 9, we found that siG3 reduced 
expression of the FR3 late protein without affecting viral DNA 

replication. To determine if the effect of siRNA to BGLF3 is due 
to its capacity to target the BGLF3 transcript, we prepared total 
RNA from samples used in figure 9. Using RT-qPCR we assessed 
die level of BGLF3 and BMRFl mRNAs in 2089 cells transfected 
with wild-type ZEBRA in the absence and presence of siG3. 
Expression of ZEBRA (Z) stimulated the BMRFl transcript by 
887-fold and the BGLF3 transcript by 41 -fold relative to cells 
transfected with empty vector (GMV). Go-transfection of siG3 had 
no effect on the level of the BMRFl transcript but markedly 
reduced the level of BGLF3 mRNA by 2.7-fold. These results 
together with data in Fig. 9 demonstrate that the observed effect of 
siG3 is due to reduction in amount of the BGLF3 transcript. 
(TIF) 

Figure S6 BGLF3 is necessary for expression of late 
genes. To investigate the effect of silencing BGLF3 on expression 

of late genes other than hfrf3 (Fig. 9A), we used RT-qPGR to 
assess the level of four late transcripts, BDLF 1 , BdRF 1 , BLLF 1 , 
and BcLFl, in the same RNA samples studied in Fig. S5. We 
found that knockdown of BGLF3 reduced the level of these late 
transcripts in 2089 cells transfected with ZEBRA plus siG3 relative 
to cells expressing ZEBRA alone by 11.2, 9.1, 10.1, and 5.9-fold 
respectively. In the same experiment (Fig. 9), the level of BGLF4 
protein was not affected by siG3. These results demonstrate that 
BGLF3 plays an essential role in regulating expression of late 
genes during the EBV lytic cycle. 
(TIF) 

Figure S7 Expression of BGLF3 and BGLF3.5 is not 
sufiicient to activate late gene expression in the absence 

of BGLF4. We examined whether compensatory expression of 
BGLF3 and BGLF3.5 in the absence of BGLF4 would rescue 

expression of the late FR3 protein. 2089 cells were transfected with 
the indicated ("xpression vectors with and without siG4. Cells were 
harvested after 48 h and expression of lytic proteins was analyzed 
by Western blot. We found that knockdown of BGLF4 abolished 
synthesis of FR3 (compare lane 2 with lane 4). Ectopic expression 
of BGLF3 and BGLF3.5 in 2089 cells with siG4 faUed to restore 
expression of the late FR3 protein (compare lane 3 with lane 7). 
This finding indicates that BGLF3 and BGLF3.5 do not substitute 
for the role of BGLF4 in regulation of late gene expression; 
BGLF4 has an independent role in activation of late gene 
expression. 
(TIF) 
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